Dragging a polymer in a viscous fluid: steady-state and transient 
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We study the conformation and dynamics of a single polymer chain that is pulled by a constant 
force applied at its one end with the other end free. Such a situation is relevant to the growing 
technology of manipulating individual macromolecules, which offers a paradigm research for probing 
far-from-equilibrium responses of long flexible biological polymers. We first analyze the Rouse model 
for the Gaussian chains for which the exact analytical results can be obtained. More realistic features 
such as the finite extensibility, the excluded volume and the hydrodynamic interactions are taken 
into account with the help of the scaling argument, which leads to various nontrivial predictions 
such as the stretching-force-dependent friction constants. We elucidate (i) generalized dynamical 
equations of state describing extension/friction laws in steady-state and (ii) the tension propagation 
laws in the transient process. We point out that the time evolutions of the dynamic friction in 
the transient process crucially depend on the experimental protocol, i.e., either constant force or 
velocity ensemble, which might be detectable in experiments using giant DNAs and chromosomes. 

PACS numbers: 82.35.Lr, 36.20.Ey, 83.50.-v 



I. INTRODUCTION 

A soft response to external forees is a generic prop- 
erty of long flexible polymers. Take a simple example 
of the static stretching, where a polymer is pulled by 
a force / [l|-ll|- Our polymer is composed of iV » 1 
segments with their size a, and its overall spatial extent 
is i? ~ aN'^ at equilibrium with the Flory exponent v. 
In good solvent condition, polymers swell, i.e., v ~ 0.6 
due to the excluded-volume interactions, where it is well 
known that a linear response domain is small, and rather 
weak force / > ksT / R suffices to enter the nonlinear 
regime {ksT is thermal energy) |l|-|3|. For stronger force 
/ > ksT/a, the inextensibility of the chain backbone in- 
troduces an additional nonlinearlity [2, Q . These aspects 
have been extensively studied both theoretically and ex- 
perimentally. 

Even more dramatic would be the dynamics, where 
polymers may exhibit intriguing nonequilibrium response 
behaviors both in a steady-state 04^] and a transient 
process piil - [l5| . This occurs when the rate of the ex- 
ternal operation exceeds the inverse of the longest relax- 
ation tinre in the system. As the relaxation times in the 
polymeric system are generally long, and their spectra 
are widely distributed, one expects that such a situation 
would be rather easily realized. The criterion for the 
onset of nonequilibrium response has the same scaling 
form as that for the significant deformation in the static 
response / > ksT/R (sec Sec. Ill Bp . 

In the present article, we study a paradigmatic case 
of the dynamical response, i.e., polymers dragged by its 
one end. When a constant pulling force is applied to 



one end of the chain at equilibrium, the chain starts to 
be deformed and settles in the steady-state conformation 
after some transient period. We analyze this "creep" ex- 
periment in a single chain level by two approaches. In 
the first, we employ a Rouse equation of motion which 
provides us with exact results. The scaling approach, 
on the other hand, allows us to unveil nonlinear effects 
associated with excluded-volume and hydrodynamic in- 
teractions while leaving unknown numerical coefficients. 
These complementary approaches offer a comprehensive 
picture of the problem. 

In Sec. nil we first analyze the steady-state, and sum- 
marize the nonlinear force-extension and force-velocity 
relations as the generalized dynamical equation of states. 
In Sec. Illli we investigate the transient process in detail 
by looking at the tension propagation dynamics, where 
a nontrivial effect of the force magnitude depending on 
the dissipation mechanism is highlighted. An additional 
insight is provided in Sec. IIVI where we formulate the 
transient process as the one-dimensional nonlinear diffu- 
sion process of local segment density towards the stretch- 
ing direction. Finally, we summarize in Sec. IVl with some 
remarks and perspectives. 



II. STEADY-STATE PROPERTIES 

We first analyze the steady-state conformation of a 
flexible polymer chain that is forced to move by a con- 
stant force applied at its one end with the other end free. 

A. Rouse model 
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Parameters characterizing our polymer chain are the 
total number A^ of bonds with equilibrium distance a 



between neighboring beads (bond length or persistence 
length). We consider the Rouse dynamics in which the 
surrounding fluid provides a local isotropic drag propor- 
tional to the local velocity, fdrag = 7r, where the friction 
coefficient 7 is proportional to rja [rj is the viscosity of sol- 
vent) and the dot represents time derivative throughout 
this paper. Since the entropic elastic forces exactly bal- 
ance the drag forces in the Rouse dynamics, the equation 
of motion is given by [l6| 



7r(n,t) 



d^r{n, t) 



-|(n,i) + f(n). 



(1) 



where the spring constant is given by fc = 3kBT/a^, 
and the random force ^ (with its mean zero) describes 
the coupling to the thermal bath and satisfies the 
usual fluctuation-dissipation relationship (^„(i)^?n(^')) = 
2^kBT5n.m5{t — t'). The constant external force applied 
at its one end s ~ Q (and along x direction) is explicitly 
described in Eq. ([T]) as 



2f5{n) 



(2) 



rather than appearing in the boundary conditions for 
r(n, t) given below. Note that the factor 2 in Eq. ^ 

assures that /^ dn f = fe^. The boundary conditions 
for r are thus those for the force-free at both ends, 
dr{n)/dn\„=o = dr{n)/dn\n=N = 0. 

Using the standard procedure based on the normal 
mode analysis, various quantities can be computed an- 
alytically [la|. The center of mass position rdt) — 
N^^ L r{n,t)dn immediately gives the force-velocity 
if-V) relation with the friction coefficient Frousc = -^7; 



^-<'°>-W 



(3) 



The polymer moves at a speed V = f /{N^) as a whole 
at steady state. The corresponding chain deformation is 
manifested in the higher modes, which leads to 
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From this, one can deduce the force-extension (f-L) re- 
lation 



L-|(r(0)-r(7V))| 



ekeT 



f 



(5) 



By comparing eq. ([S]) with the equilibrium coil size R = 
aN^^^, one indeed find the threshold force / = 6kBT/R 
above which the deformation becomes apparent. One 
can also obtain the velocity-extension (V-L) and f-V-L 
relations by eliminating f or N from eq. ([3]), (O as 
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(6) 



The dragging force builds up along the chain from the 
downstream end; the internal tension at length n in the 
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FIG. 1: Schematic representations of a dragged Rouse chain 
by a constant force /. The bottom picture corresponds to the 
case with high force f > fc where there is a stem of length 
ma in the front region. 



polymer chain is along the x direction, T(7i) = f{n)ex, 
where f{n) = k{dr/dn) = /(I — n/N). 

Second moment: Second moment quantities allow us 
to gain more information on how the stretching alters 
entropic properties in our Gaussian random coil confor- 
mations. A straightforward calculation yields 

(|r(?i) — r(?7i)p) = a{n — m) 
t2„2 / /• \ 2 r 
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(7) 



The first diffusive term describes entropic randomization 
effects while the second term arises from the stretching 
of the chain. Comparing these two contributions leads 
to the characteristic size A = n — m (as function of n) 
separating entropy-dominated and stretching-dominated 
regimes. Assuming a sufficiently long chain iV 3> 1, we 
obtain, up to the first order of IS./N , {af{n)/kBT)^A ~ 
1. where f{n) = /(I — n/N) is the tension at n given 
above. Introducing a n-dependent blob size as £,{n) ~ 
aA^/^, we find the relation 



f{n)an 
SksT 



1. 



(8) 



This has a clear and important physical meaning (see 
Fig. [1] (top)); Sitting on the position of the bead n, and 
look at the chain conformation around it. The tensile 
effect will be a weak perturbation and the chain confor- 
mation is essentially the same as the trajectory of a ran- 
dom walk at the scale smaller than ^(n) ~ 3kBT/f{n). 
At larger length scales, however, the chain is highly de- 
formed to the stretched conformations. An asymmetric 
way of the driving, i.e., pulling one end, leads to the 



nonuniform tension profile and the size-dependent blob 
size £,{n). 

Effect of the finite chain extensibility: As we have seen, 
the tension is highest at the pulling site. From eq. (|H), 
the distance between the dragged bead (n = 0) and the 
next bead (n = 1) is calculated as 



(|ri-ro|) 
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(9) 



From this, we notice that when the force reaches the 
threshold 



/c = 



SfcsT 



(10) 



the bond reaches its maximum extension (|ri — tqI) = a. 
For stronger forces, the Rouse model which maintains 
the bond by a harmonic potential is no longer applica- 
ble. Below, we will discuss this effect of the finite chain 
extensibility in a simplified manner. 

Consider that the chain is completely extended up to 
the m-th bead counting from the dragging site, forming a 
stem of the length ma (see Fig. [T] (bottom)). This means 
that the tension at the m-th bead is equal to fc- For the 
remaining rear part with N ~ m segments, the results 
for the Rouse model obtained above may be applied, so 
that using Eqs. ([5]) and (|T0|) . we obtain its extension ~ 
a^ {N —m) fc / QksT = {N~m)a/2. Since the extension of 
the stem is ma, the total end-to-end distance of this stem- 
flower shape is L = (TV -I- m)a/2. The number m can be 
determined from the total force balance / = N'jV ~ fc+ 
mjV, leading to m = {f-fc)/{lV) = N{l-3kBT/fa). 
The force-extension relation in this strong force regime 
is 
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(11) 



It is instructive to compare these results eq. (O 
and pT|) with the force-extension relation for the free- 
jointed chain pulled by its both ends, i.e., a static stretch- 
ing (Fig. [2]). In this case, the tension thus the confor- 
mation is uniform along the chain, and the exact force- 
extension in the full force range is available [2]: 
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(./ < ksT/a) 
if > keT/a) 



(12) 



Recall that a freely-jointed chain and a Rouse chain arc 
similar in the sense both have no excluded volume effect, 
but the former is now in the static situation pulled by 
both ends while the latter is in the dynamical steady state 
with strong inhomogeneity along the chain. Nonetheless, 
we see that there is no distinction between these two in 
the scaling level. Difference, however, clearly exists in 
the exact prcf actors, i.e., the one-site dragging results in 
less extension in the same force as is expected. 
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FIG. 2: Force-extension relations of a dragged Rouse chain 
pulled by its one end and a free-jointed chain pulled by its 
two ends. 



We have adopted the simplest prescription to analyze 
the finite chain extensibility effect at high force f > fc- 
In general, however, more sophisticated approaches re- 
quire the specification of the fiexibility mechanism in the 
segment scale. This can be, for instance, a freely-jointed 
model or a work-like-chain model which is suitable for 
the chain with a uniform bending elasticity. The present 
simple prescription, however, suffices for the investiga- 
tion of scaling properties. In the following scaling ap- 
proach, we shall see that both cases (/ < fc and / > fc) 
can be treated in the same formalism, but with different 
characteristic scaling exponents a and /3 defined below 
(eqs. ^-^). 



B. Scaling approach 

The Rouse model is applicable to polymers in 
melts [l6|. But in many of other situations, the excluded- 
volume interactions and/or the solvent mediated hydro- 
dynamic interactions should be taken into account. The 
former gives rise to long range correlation along the chain 
and the latter makes the friction coefficient conformation 
dependent. The resultant nonlinear effects can be treated 
within the framework of the scaling theory. 

Onset of the strong deformation: Assume that a 
pulling force is weak enough so that the dragged chain 
takes an equilibrium conformation R ~ aN'^ with the 
Flory exponent i>. We denote the friction coefficient of 
the chain as F = j{R/a)^~'^ with the so-called dynamic 
exponent z [l|. A longest relaxation time of the poly- 
mer coil is then r^q ~ TR^/ksT ~ To{R/a)^, where 
To = ja'^/kBT ~ Tja^/kBT is a segment scale micro- 
scopic time. Comparing this with the typical velocity 
gradient 7 ~ V/R. ~ {R/ay~^ f/ia-f), one finds a char- 
acteristic force / ~ kBT/R above which the chain ex- 
hibits substantial deformations. Note that this coincides 



with the static criterion given by comparing sizes of the 
Pincus blob ksT/ f and the equilibrium coil R [iH3- For 
weaker force, the chain is assumed to be in near equilib- 
rium where ordinary linear response applies. 

Steady-state conformation: The conformation of a 
tethered chain submitted to a solvent flow of velocity 
V was studied by Brochard-Wyart using scaling ap- 
proach [5|, |6| . Here we present its generalized form using 
two critical exponents v and z. Let us assume that the 
force is moderately strong, i.e., ksT/R < f < ksT/a ~ 
fc- As we have already seen in the analysis of Rouse 
model, for a chain dragged in a viscous fluid, there is 
a blob size £,{x) below which the effect of the pulling is 
insignificant, and an equilibrium formula for the Flory 
relation ^{x) ~ ag{xY can be applied, where g repre- 
sents the number of segments in the blob. Note that the 
one-dimensional coordinate x is taken along the long axis 
of a stretched polymer chain (whose shape is assumed to 
have an axial symmetry) , and not the arc- length param- 
eter as in the previous section. We assume that the same 
scaling behavior holds for the local friction coefficient 



r(0 



(13) 



Hydrodynamic interaction among beads within a blob 
can be taken into account by setting z = 3 (non-draining 
case), leading to the Stokes formula F ~ 77^. On the other 
hand, setting z = (1 + 1v)lv amounts to assume that the 
solvent is just immobile (free-draining) . We then recover 
the Rouse-type friction law F(^) ~ r]ag. The overall de- 
formed conformation can be pictured as a sequence of 
blobs of size ^ which are hydrodynamically decoupled. 
Hence, the hydrodynamic friction force acting at the po- 
sition X can be evaluated as an integral of the contribu- 
tion which builds up from the downstream free end; 



/drag (a;) ~ V 



dx — 



(14) 



Applying a local force balance equation ^ 
ksT/ f drag 13|, we obtain the conformation profile as 



a^) 



TqVx 



ll(1-z) 



(15) 



By equating the friction force at the pulled end to the 
applied force (total force balance) , we find a steady-state 
f-V-L relation, i.e., a generalized version of eq. ^ in the 
form 



f^ksT/roVLY^^-'^ 



(16) 



To obtain a generalized V-L relation, we integrate the 
number of segments from the downstream free end to the 
pulled site, and equate it with the total segment number, 
i.e., mass conservation law: TV = /^ {g/£)dx, which leads 
to 



a 



(v-l)/(l+v-zv) 



(17) 



It is convenient to summarize eq. p6|) and eq. p7)) as the 
following relations which we call the dynamical equation 
of state; 
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(18) 



(19) 



(20) 



where defining two exponents 



z - 2, /3 = (1 - v)/v. 



a = z - 'Z, p = (L - v}/v. (21) 

turns out to be useful, but it should be kept in mind 
that these are merely transcriptions of the original ba- 
sic dynamic and static exponents z and v. Note that 
all of these relations reduce to those of Rouse model 
for an ideal statistics v — 1/2 in the free-draining limit 
z = (1 -|- 2i/)/y = 4. However, the excluded-volume and 
hydrodynamic interactions alter the exponents. In par- 
ticular, one can see the nonlinear friction law of the hy- 
drodynamic origin from eq. (J19p and the nonlinear ex- 
tension law of the excluded-volume origin from eq. (j20p , 
i.e., the linear f-V or f-L relation holds only for the 
free-draining condition or for the chain with v = 1/2, re- 
spectively. Note also that the scaling f-L relation is the 
same as that for a well-known situation of the static chain 
stretching, i.e., pulling force applied to both its ends [l|- 
y|. Therefore, the difference in the force-extension re- 
lation between dynamic and static stretching is mani- 
fested only in the numerical coefficient (compare eq. ([5]) 
and (gni) with v == 1/2). At the free end of the chain, 
the largest blob of size ^f,oc there experiences the hydro- 
dynamic dragging force acting on itself only. From the 
Pincus relation fesT/^froc — 7^(^frcc/a)^~^, we find 
/_ T^x 1/(1-^) 

6ree - « ^ (22) 



Effect of the finite chain extensibility: So far we have 
dealt with the "trumpet" regime f < fc- For stronger 
force, we have to take account of the finite chain exten- 
sibility, as we have seen in the analysis of Rouse model 
(Sec. Ill Ap . where the dragged chain takes a so-called 
"stem- flower" conformation [6|. Our formalism is appli- 
cable even to such a situation, i.e., a set of basic equations 
eqs. (|18p - (|^D|) remains intact with the assignment of ap- 
propriate exponents a = 1 and /3 = 0. Indeed, /3 = 
in eq. (PD|) and a = 1 in eq. p^ indicate L ^ N and 
/ ~ A^T^, respectively, i.e., full stretching on the scal- 
ing level Q (see eq. pT|) for Rouse model). The size of 
the largest blob ^ftoc is still given by eq. p2|) until the 
point ^frcc > a <^ f < NksT/a, where use has been 
made of eq. (fT9|) with a — I and /3 = 0. For stronger 
force / > NksT/a, an entropic coiling at the rear end 
is negligible so that the dragged chain takes almost fully 
stretched conformation. 



III. TRANSIENT DYNAMICS 

Now we turn our attention to the transient dynam- 
ics after the sudden action of a puhing force at one end. 
Here we wih observe that there are a variety of tension- 
propagation dynamics depending on the types of hydro- 
dynamics and/or excluded volume interactions. 



A. Rouse model 

To set the stage, we again start with the analysis of the 
Rouse model dynamics. The basic equations are already 
given in Sec. Ill Al Wc only need to modify cq. ^. The 
dragging force is switched on at t = 0, then, 



/„ = 2fuit)Sin)e, 



(23) 



where u{t) = 1 for t > and u{t) = for f < 0. The 
resulting linear equations can be solved exactly again us- 
ing the normal mode analysis. We are particularly in- 
terested in the apparent pulling velocity of the polymer 
chain (i.e., the velocity of the pulling end), V{t) ~ (ro), 
which is found to be 



V{t) = ^[1 + G(i)]e„ 



where the propagator G{t) is 



G{t) 



2^exp 
p=i 



7^ 



pH 
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(24) 



(25) 



In the last near-equality, we have replaced the sum- 
mation with the integral, and the Rouse time tr = 
-fa^N'^/iSnkBT) has been introduced ^. The final re- 
sult 



Vit) 






-1/2 



(26) 



indicates that in early stage t ^ tr, the velocity de- 
creases as V{t) ~ f/{Nj){t/Ti{)^^/^ex, and eventu- 
ally reaches the steady-state with the velocity V{t) = 
f /{N^)ex after the characteristic time tr. Its physical 
interpretation is as the following. At i < tr, only a part 
of segments M{t) < N close to the dragged site can re- 
spond to the pulling force, while other N—M{t) segments 
in the rear are essentially unaffected yet, thus taking an 
unperturbed conformation at rest. Then, we have the 
force balance equation, ^M{t)V{t) ~ /, from which we 
obtain 



Mit) 



1/2 



(27) 



The tension applied at the dragged site diffusively prop- 
agates along the chain backbone. 




-R (= -aN') 



FIG. 3: Schematic representations of the transient dynamics 
in stretching process. 



B. Scaling approach 

Imagine that one end segment located at the origin is 
started to be pulled to the x direction at i = 0. Be- 
fore the pulling force / is applied (i < 0), the chain at 
rest takes the equilibrium conformation. As already dis- 
cussed in Sec. IIIBl the whole chain would follow such 
a force with a velocity V = //F by retaining equilib- 
rium conformations, if the force is below the thresh- 
old / ~ ksT/R. However, for stronger forces, only 
the fraction of the chain (with the number of segment 
.90 ^ {(.a/af^" ^ {fa/kBT)-^'") close to the driving site 
can immediately respond and set in motion, while keep- 
ing the remaining rear part at rest yet. Such a transient 
response can be described by what we call the two-phase 
formalism originally introduced in the problem of driven 
translocation across a nanopore |12l . Il3| , in which one an- 
alyzes the dynamics of the front separating the respond- 
ing/unresponding domains, in other words, the dynamics 
of the tension propagation along the backbone. 

At time i, the positions of the front and the dragged 
segment are x = —X{t) and x = ^i), respectively 
(see Fig. [3]). The responding part composing M[t) seg- 
ments in the range x £ [—X{t), l{t)] assumes steady-state 
conformations moving with the representative velocity 
~ V{t), for which the result obtained in Sec. Ill Bl can 
be applied with the suitable redefinition of the variables, 
i.e., L ^ 1 + X, N ^ M. From cqs. ^ and dH]), we 
obtain 
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The moving distance of the driving site is 



m 



<\sV{s) 



(30) 



Recahing that the chain is essentially unaffected behind 
the front, one can notice that the front position X{t) and 
the segment label M{t) (counting from the pulled end) 
at the front satisfies the relation 



X{t) ~ aM". 



(31) 



Equations ((28l) - ((3T|) constitute a closed set of the rela- 
tions of V{t), M{t), R{t). The solution can be written in 
the form 



QiV{t)) - QiVo) = - 

To 

where we have defined the function 
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with Kis being the steady state velocity given by eq. ([T9| . 
The last near-equality in eq. p3p holds for the situation 
of our interest / :^ fc^T/i?. 

Except for the very initial period, eq. p3p can be ap- 
proximated as Q{V{t)) ~ i/ro which can be arranged 
into 
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(34) 



where the characteristic time of the transient period can 
be deduced from the condition Ttr/ro ~ Q{V{Tt^)) = 
Q{Vss) so that 



Ttr ^ TqN' 



f fa 



\kBT 



(35) 



with the exponent 



2/3-Q = 



{2/iy) - z {knT/R <f< kBT/a) 



-1 



(/ > ksT/a) 



(36) 



Note that in the moderately driven regime {ksT/ R < 
f < kBT/a), eq. ([55]) can be rewritten as Ttr — 
TeqifR/ksT)'^- At the threshold force / ~ ksT/R, we 
see the crossover between the non-equilibrium and equi- 
librium dynamics, Ttr — Toq, where r^q = To{R/a)^ = 
tqN'^^ is the longest relaxation time of the equilibrium 
coil. 



Combining cq. (|34)) with eq. ((29)) . we obtain the time 
*> '' evolution of the propagation front as 
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fa 
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(37) 



We thus find that the tension always propagates diffu- 
sively along the chain backbone with the force dependent 
amplitude. The exponent uj determines the force depen- 
dence in the transient process (eqs. ([55]) and ([57)) ). For 
Rouse chain in the trumpet regime (v = 1/2, z = 4), we 
have oj = 0, in agreement with the exact result obtained 
in Sec. IIII Al With the excluded- volume effect included 
{v > 1/2) in the free-draining case, we have the force 
amplitude exponent uj = (1/j^) — 2 < 0. Surprisingly, 
this exponent reverses the sign in the non-draining case; 
ijj ~ (2/t^) — 3 > irrespective of the excluded-volume 
effect. One can obtain a clue on this somewhat counter- 
intuitive result by rewriting eq. ()35p as 
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N~\fa/kBT) 
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(38) 



where we have utilized the dynamical equation of states, 
i.e., eqs. (fT9)) and (|20)) . While both the steady-state ex- 
tension Lss and the velocity Vss increases with /, their 
relative rate generally depends on the dissipation mech- 
anism and the excluded-volume effect. 

Under the stronger force (/ > ksT/a), the chain is 
highly stretched so that the excluded-volume effect be- 
comes irrelevant, and we would expect the exponent 
cj = — 1. It may be worth noting that the initial ve- 
locity Vb is given from the condition r[^o]Vb = /; with 
the initial blob size ^o — ^sT// 
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This has the same scaling structure for the downstream 
blob size ^froo in the steady state (eq. (|^ ). Just as ^froo 
saturates to the segment size at large force / > NkBT/a, 
the initial blob size reaches the lower bound S^q :^ a at 
the threshold force / ~ kBT/a, thus, the initial velocity 
is modified as Vb — {a/To){fa/kBT) for larger forces. 



IV. NONLINEAR DIFFUSION EQUATION 

The above scaling predictions based on the two-phase 
formalism can be corroborated by a semi-quantitative 
argument that considers a viscous dynamics of the lo- 
cal segment density in a single polymer chain. We will 
here outline this approach and derive a nonlinear diffu- 
sion equation that can describe the spatial profile of chain 
deformations in the transient process. Analytic solutions 
are available to our nonlinear diffusion equation, from 
which a number of the scaling predictions are confirmed. 

Take one dimensional coordinate x along a pulling di- 
rection. The cross section of the polymer at x is ^{x,t). 



the velocity field of the segments averaged over this cross 
section is mainly along x direction, which we denote as 
v{x,t). The mass conservation averaged over the cross 
section at x can be written as 



|«^) + ^(.f'.) 



0, 



(40) 



The segment density p satisfies p^^ ~ a^g, where g is the 
number of segments within a blob of size ^, and could be 
given due to the Flory relation as f ~ o-g'^, from which 
we obtain p^^ ~ a?[S^/aY . The local force balance along 
x-direction may be written as 



dT 
dx 



r 

T 



-V = 0, 



(41) 



where T(x, t) is an internal tension in the polymer chain 
(averaged over the cross section at a;), and T is given 
in eq. ([T^. The Pincus rule suggests T ~ kBT/£_. 
The boundary condition at the pulling end, x = l{t), is 
T{l,t) = /, or using T ~ ksT/^, it gives on the scaling 
level the boundary condition for ^ as 



ai,t) 



f ■ 



From Eqs. (|40| and p4|) . we obtain 
w(x, t) 



ksT f a\^ d fa 



1 



dx VC 



(42) 



(43) 



Substituting this into eq. (|40p and introducing the line 
segment density (f) — ag/(, ~ (^/a)*^, we arrive at the 
nonlinear diffusion equation 



dt 



d_ 

dx 



h-^^/p: 



dx 



(44) 



where the diffusion constant k includes unknown numer- 
ical factors which ensure the equality in Eq. pi)) . i.e., 
K — (const) X P^^ksT/j. 

Solution and scaling behavior- Equation (|44p or its 
analogue has recently been proposed to describe the 
polymer detachment kinetics from absorbing surface [l5| 
or the decompression/unfolding dynamics of an initially 
compacted polymer [l^. More traditionally, this type 
of nonlinear diffusion equation is often encountered in 
the ecology, i.e., insect and animal dispersal and inva- 
sion [l3l. It also describes other various phenomena, 
sometimes called as the fast {—a/ (3 < 0) or the slow 
{—a//3 > 0) diffusion process, the relevance of which 
ranges from plasma physics, kinetic theory of gases to 
the fluid transport in porous medium. 

The scaling relation for the dynamics of pulled end 
simply follows from the self-similar property of the sys- 
tem. Let 1{t) and (J)o{1{t),t) be the position of the 
pulled end at time r and the segment line density there. 
The self-similarity implies that eq. (j44p is invariant un- 
der the scale transformation t ^- rt, x ^>- l(t)x and 



(p -^ <t>o4>- This leads to the relation P ~ tk^q " 
This, together with the boundary condition (eq dH])) 
00 = [fo-/kBT)~^ , leads to the time evolution of the 
pulled end l{T)la ~ {r/ay/'^ifa/kBT)"/'^, which is pre- 
cisely in accord with eq. (p4|) . 

To obtain a solution, it is convenient to introducing the 
variable h = (^/a)~°. Then eq. (|44p can be rewritten as 



dh 
'dt 



Pk 
a 



dh 
dx 



K,h 



d^h 
dx"^ ' 



(45) 



Analytic forms of self-similar solutions for eq. pS)) are 
known for the case a//3 < 2 [l7|- The special solution 
that can satisfy the boundary condition (|^^ is h{x^ t) = 
A{t) + B{t){x/af., where 



Ao (1 + cjot)"/(2/J-") (46) 

Bo{\+ujQtr\ (47) 

1 
The polymer shape in the transient is thus found to be 

l/a 



A{t) 
Bit) 

with LOo = 2BoKa~^[{2/3/a) - 1] 



ax,t) = ao,t) 



1 



A2(i) 



(48) 



where ^(0,i)/a = [A{t)]-^^°', and X{t)/a = 
[A(i)/B(t)]i/2 - [1 + ujotf^'-^^-^K From this and 
eq. (|43p . the velocity field is obtained as 



v{x,t) c^ B{t) — . 
To 

This allows us to write the profile as 

^(x,t)/a^[B(i)(x/a)2]-i/"-^''°''^^''^- 



(49) 



-l/a 



(50) 



which is valid in the region a; ^ A. This coincides with 
the form of steady state profile given in eq. ([T5]). 

The boundary condition, eq. (|^ . requires A{t) + 
B{t){l{t)/af ~ {fa/ksT)"', from which we obtain the 
time evolution of the position of the pulling end l{t). In 
particular, when (fa/kBT)"' ^ A{t), that is, when 



//i?\('^-"' 
^/- « [l^) 



(51) 



where R ~ ^(0,0) is a typical equilibrium size of the 
polymer chain, we can find the scaling law for l{t) given 
by 



lit) 



fa 
keT 



a/2 



1/2 



(52) 



which again accords with eq. ([M)) . This scaling behav- 
ior is expected to be observed during the time region 
given by 1 « t/ro < {fR/kBT)^^'^-"^ ^ nr/req. There 
will be a substantial time range satisfying this, given the 
nonequilibrium regime fR/kBT ^ 1 under considera- 
tion. 



As is evident from the time dependency of A, the 
method of present analysis yields a physically meaningful 
solution when a//3 < 2. This condition is always satisfied 
for the non-draining dynamics {z = 3), but not for the 
free-draining dynamics {z = [2v + \)lv). In this regard, 
it is interesting to note a recent study on the polymer de- 
tachment kinetics, which demonstrates Z(i) ~ i^/^ scaling 
for the pulled end even when a//3 > 2 [15| . Although the 
adopted linearizing approximation of eq. (j44p in ref. [l5| 
limits the range of their analysis to the weak force case 
only, it provides a further support for the scaling anal- 
ysis in Sec. IIIIBI We also note that the point 2/3 = a 
corresponds to the Rouse chain, i.e., z = 4, i^ = 1/2 and 
coincides with the condition w = at which the force 
exponent reverses its sign (cq. ([36| ). The Rouse chain is 
exactly analyzed using the Gaussian model in Sees. Ill Al 
and IIII Al and the results arc already confirmed to agree 
with scaling predictions in Sec. IIIIBI 

Lastly, we point out that the results obtained in this 
section shed some light on the two-phase formalism de- 
veloped in Sec. IIIIBI The linear velocity profile eq. (|49|) 
leads to the average velocity of the moving domain 

Vave = SxtZ"" 'l>{x)dx/{xhead " Xfail) = {vhead " Vtail)/'^, 

where subscripts "head" and "tail" indicate the begin- 
ning and the end points of the moving domain. Given 
Vtaii ^ Vhead, it is natural to set the representative veloc- 
ity as T^ ~ Vhead, which corresponds to our equation ([30|) 
in Sec. IIIIBI The spatial profile of the moving domain in 
the transient process is sharper than that in the steady- 
state, i.e., compare eqs. ([50]) and ([T5|) because of the ve- 
locity gradient in the moving domain in the former situ- 
ation (eq. dH])). 



V. SUMMARY AND PERSPECTIVES 

We have presented a comprehensive analysis for the 
dynamics of a dragged polymer, a paradigm for the 
nonequilibrium response of polymer chains. Both the 
steady-state and transient properties have been clarified 
by means of the exact solution of Rouse model and the 
scaling approach. The latter allows us to establish the 
nonlinear extension/friction laws in steady-state and the 
tension propagation law during the transient process in 
general form, which involves the Rouse model result in a 
particular case. 

The present set-up of dragging chain could be real- 
ized by manipulation techniques using optical/magnetic 
tweezers or atomic force microscopy (AFM). Here, there 
are two basic protocols for the manipulation, (i) A con- 
stant force is applied to the end segment, and (ii) the 
end segment is dragged in a constant velocity. This lat- 
ter case is realized by fixing the end segment in space 
and then moving the stage in a constant speed. During 
the transient period the friction constant increases with 
time, and our analysis indicates that the time evolutions 
of the dynamic friction are generally different in these two 
ensembles. In the constant force ensemble, the friction 



constant grows as 



r/(t) 



v{t) 



7 



fa 



^1 / l\ "^z.* 



T-Q 



(53) 



with exponents (5/ = 1 — a/2 and 5f,t = 1/2. On the 
other hand, in the constant velocity ensemble, then, the 
evolution of the friction constant is 



Tv{t) 



fit) frovywt 



Sv,t 



(54) 



with Sv = (2/a) — 1 and 6v,t = 1/a. 

Transcribing back to the original dynamical exponent 
in the trumpet regime k^T/R < f < k^T/a, we have 
<5/ = 2 - (z/2), 5v = (4 - z)/{z - 2) and 6v,t ^ l/(^ - 2). 
For a chain with ideal statistics ly = 1/2 in the free- 
draining limit z = [1 + 1v)lv = 4, there is no distinc- 
tion between these two expressions, r/(t) = Vyify — 
^{t/ToY^^, which coincides with the result obtained for 
a Rouse model in cq. (j27[) . In other cases, hydrodynamic 
and/or excluded-volume interactions give rise to distinct 
dynamical behaviors for each of two ensembles character- 
ized by nonzero Sf, Sy as well as the difference in Sf^t and 
6vt- In particular, (i) in the non-draining approximation 
(z'= 3), Sf = 1/2, 6v = 1 and 6f.t = 1/2, Sv,t = 1; 
(ii) in the free-draining approximation (z — {I + 2v)/v)^ 
6f = I- (1/2j/), 5v = 2v -I and 5f,t == 1/2, 5v,t = v- 
In both cases, the time evolution arc different depend- 
ing on the ensemble. For stronger force k^T ja < f, the 
chain is highly stretched, thus the result becomes insen- 
sitive to excluded- volumc/hydrodynamic effects. Substi- 
tuting a = 1 and /3 = Q, we find 5f = Sf^t = 1/2 and 
5v = 6v,t = 1- In usual experiments in aqueous media, 
where the non-draining approximation is expected to be 
valid (neglecting a logarithmic correction), we thus pre- 
dict Sf = Sf,t = 1/2 and Sy — Sy.t — 1 in the entire 
force range / > k-QT/ R. We have to keep in mind, how- 
ever, that a very long chain would be required to obtain 
clear signals by avoiding the effect due to a micron-sized 
trapping bead in tweezers experiment and the cantilever 
in AFM. Giant DNAs and chromosomes (appropriately 
treated if necessary) might be good candidates to test 
our predictions. 

Closely related to the present analysis is the dynam- 
ics of polymer translocation through a narrow pore, on 
which there have been numerous recent researches [20[. 
It has been proposed and verified that the dynamics 
of driven translocation is described as the sequential 
nonequilibrium process associated with the tension prop- 
agation along the chain [l3, [IjI akin to the transient 
process of the dragged chain. Only distinction lies in 
the way how the driving force is exerted, i.e., unlike the 
end-pulling for the dragged chain (eq. (Q])), the force is 
acting only on the segment inside the pore which is fixed 
in space. Since the polymer exhibits a biased Brown- 
ian motion, the segment under the action of the driving 
force changes with time. This makes a trouble in the im- 
plementation of the force term in Rouse model, and also 



causes some subtlety in the description of the moving do- 
main in the scahng approach [2l], [g^]- Nonetheless, the 
basic physics can be well captured by the two-phase for- 
malism, which can be developed along the same line with 
the present problem. We expect that a list of related ex- 
amples would build up by adding phenomena in biology 
and polymer science, where the notion of the noncquilib- 
rium response plays an important role. 
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